Polar distortions in hydrogen bonded organic ferroelectrics 



Alessandro Stroppa,^ Domenico Di Sante,^ Sachio Horiuchi,^ 
Yoshinori Tokura,^ David Vanderbilt,^ and Silvia Picozzi^ 

^CNR-SPIN, L'Aquila, Ital^ 
^University of L'Aquila, Physics Department, Via Vetoio, L'Aquila 
^National Institute of Advanced Industrial Science and Technology (AIST) Tsukuha, Iharaki, 305-8562, Japan 

and 

CREST, JST, Tsukuha, 305-8562, Japan. 
^Department of Applied Physics, The University of Tokyo Hongo Bunkyo-ku, Tokyo, 113-8656, Japan 

and 

Correlated Electron Research Croup, RIKEN Advanced Science Institute, Wako 351-0198, Japan 
^Department of Physics and Astronomy, Rutgers University, 
136 Frelinghuysen Road, Piscataway, New Jersey 08854-8019, USA 
^ CNR- SPIN, L'Aquila, Italy 
(Dated: January 20, 2013) 

Although ferroelectric compounds containing hydrogen bonds were among the first to be discov- 
ered, organic ferroelectrics are relatively rare. The discovery of high polarization at room tempera- 
ture in croconic acid [Nature 463, 789 (2010)] has led to a renewed interest in organic ferroelectrics. 
We present an ab-initio study of two ferroelectric organic molecular crystals, l-cyclobutene-1,2- 
dicarboxylic acid (CBDC) and 2-phenylmalondialdehyde (PhMDA). By using a distortion- mode 
analysis we shed light on the microscopic mechanisms contributing to the polarization, which we 
find to be as large as 14.3 and 7.0/iC/cm^ for CBDC and PhMDA respectively. These results sug- 
gest that it may be fruitful to search among known but poorly characterized organic compounds for 
organic ferroelectrics with enhanced polar properties suitable for device applications. 



The property of polarization switchable by an applied 
external electric field, i.e., ferroelectricity, is the basis of a 
wide range of device applications [l^. The first known fer- 
roelectric material, discovered in 1920, was sodium potas- 
sium tartrate tetrahydrate (NaKC4H406-H2 0), better 
known as Rochelle salt [2]. Unfortunately it is also one 
of the most complicated ferroelectrics known to date, 
and research in this field soon focused on simpler fer- 
roelectrics, such as phosphates and arsenates, discovered 
later 0,0]. A prototypical example is potassium dihydro- 
gen phosphate, KH2PO4, better known as KDP[5]. The 
latter contains hydrogen bonds for which different possi- 
ble arrangements of the hydrogens can result in different 
orientations of the dipolar units. After the discovery of 
ferroelectricity in BaTiOa, with polarizations as large as 
27/iC/cm^ in the tetragonal phase, [6] attention switched 
to this new class of perovskite oxygen octahedral fer- 
roelectrics made up from basic BOe building blocks, of 
which BaTiOa was the forerunner. These perovskite and 
related materials are by far the most investigated class of 
ferroelectrics, and the most important for current device 
applications. Nevertheless, substantial efforts are now 
being made in order to find ferroelectric materials that 
are potentially cheaper, more soluble, less toxic, lighter, 
or more flexible [t*, Very recently it was discovered 
that organic crystals of croconic acid, H2C5O5, exhibit 
ferroelectricitl y w ith a large spontaneous polarization of 
21 /iC/cm^ [9|,1Iq|. Croconic acid consists of polar stacks 
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of sheets of hydrogen-bonded molecules [11]. Upon ap- 
plication of an electric field, protons associated with one 
molecule cooperatively shift to a hydrogen-bonded neigh- 
bor, switching the molecular dipole and giving rise to a 
large polarization [9]. It seems likely that this physical 
mechanism may also occur in many other organic ma- 
terials whose ferroelectric properties have yet to be dis- 
covered and characterized. Ab-initio calculations are a 
suitable tool in this case, as they can not only identify 
materials with large polarization, but also shed light on 
the mechanisms behind the polarization itself. 

In this work, we focus on l-cyclobutene-1,2- 
dicarboxylic acid (CBDC, C6H604)[il and 2- 
phenylmalondialdehyde (PhMDA, CQH802)lil. Initial 
x-ray diffraction studies of CBDC[l2] and PhMDA jlsj 
were performed only at room temperature for CBDC 
and at T=111K for PhMDA. It was only recently that 
new measurements were extended to low temperature 
for CBDC and to room temperature for PhMDA. 
showing that both compounds are ferroelectric with 
Curie temperatures above 400 and 363 K respectively. 
The observed spontaneous polarization (Px,Py,Pz) is 
as large as (0.4,0,2.8) and (0,0,9) /iC/cm^, for CBDC 
and PhMDA respectively [14]. The main purpose of 
the present study is to describe the ferroelectricity 
in these compounds based on ab-initio calculations. 
We have found that for both materials, the main 
contribution to the polarization is a collective proton 
transfer between molecular units, hereafter denoted as 
an "inter- molecular" hydrogen shift. This is basically the 
same main mechanism found in croconic acid. However, 
importantly, we have also introduced a partial mode 
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analysis of the contributions to the ferroelectricity, and 
found that correlated intramolecular distortions also 
contribute significantly to the polarization. 

First-principles density functional theory calculations 
were performed usin g th e Vienna-Ab-Initio Simulation 
package (VASP)p!5l, Tl6|. The Kohn-Sham equations 
were solved using the projector augmented wave pseudo- 
potentials and the PBE generalized gradient density ap- 
proximation to the exchange-correlation potential |17|. 
We used a plane- wave cut-off of 400 eV and k-point 
meshes of (6,2,4) and (4,2,6) for CBDC and PhMDA re- 
spectively. The Berry-phase method was employed to 
evaluate the crystalline polarization [l8[ [l9|. All atomic 
positions were optimized until the forces were below 
0.01 eV/A. Test calculations were performed to estimate 
the effect of (1) effects beyond the local density ap- 
proximation by using the Heyd-Scuseria-Ernzerhof hy- 
brid functional as implemented in VASP,[20-22] and (2) 
van der Waals corrections [23, 24] as proposed within a 
density-functional framework. In both cases the changes 
in the magnitude of polarization were found to be less 
than a few percent, confirming the basic physics ex- 
plained in the present work. 

In Fig. [T] we show a ball-and-stick model of CBDC 
(top) and PhMDA (bottom). In (a) and (c) we show 
the centric structures, and in (b) and (d) the polar ones. 
The CBDC molecular unit is formed by a main planar 
four-membered ring of carbon atoms similar to the pla- 
nar cyclobutene molecule, and two carboxyl carbons, i.e., 
with formula — C=0— O— H, where and ' — ' refer to 
double and single bonds. The molecular units shown in 
Fig. [1] are part of an infinite chain of molecules related 
by a glide plane and linked by intermolecular hydrogen 
bonds. The PhMDA molecular unit is formed by a planar 
phenyl group, i.e., six carbon atoms arranged in a planar 
ring, each of which is bonded to an H atom. The phenyl 
group is linked to a linear hydrogen-bonded chain of j3- 
diketone enol moieties, i.e., 0(=)C— C=C— O— H. Hydro- 
gen bonds Hnk the molecules together along the [102] and 
[102] directions into infinite chains. 

We found it useful here to introduce a pseudosymme- 
try analysis, in which a given low-symmetry (ferroelec- 
tric) structure is represented in terms of a symmetry- 
lowering Landau-type structural phase transition from a 
high-symmetry (paraelectric) parent structure, based on 
finding a supergroup of the given space group [251 l26j . 
CBDC crystallizes in the monoclinic space group Cc and 
its pseudosymmetric centric structure has space group 
symmetry (72/c; the two structures are related by a max- 
imum atomic distortion of 0.43 A. Analogously, PhMDA 
crystallizes in the Pna21 and its pseudosymmetric centric 
structure is Pbcn^ with maximum atomic displacements 
of 0.25 A. In order to gain insight into the ferroelectric- 
ity, we compare the relaxed centric and polar structures 
shown in Fig.^a-b) for CBDC and in Fig.^c-d) for Ph- 
MDA. For CBDC we have two types of hydrogen bonds 
("intramolecular" and "intermolecular") which, in the 
polar state, shift toward the molecular units on the right. 



FIG. 1: (Color online) Ball-and-stick models of CBDC (top) 
and PhMDA (bottom), (a) and (c) refer to centric and (b) and 
(d) to polar structures. Dashed guiding-eye lines refer to the 
position of relevant hydrogens contributing to polarization; 
arrows in (b) and (d) indicate important polar distortions. 



as shown by the short arrows in Fig. [IJb). There is an- 
other cooperative atomic distortion, shown in Fig. [Hb) 
by the curved arrow, hereafter referred to "molecular 
buckling." For PhMDA, the hydrogen sitting between 
two molecular units in the polar structure shifts towards 
one of its neighboring units, as indicated by short arrows 
in Fig.[T]^d). Two other relevant atomic distortions come 
into play, as shown by the curved arrows, both tending 
to deform the molecular units. One acts on the planar 
phenyl group, while the other distorts the /3-diketone enol 
moieties. For both compounds, then, we find that three 
different types of distortion can contribute importantly 
to the ferroelectric polarization. 

In Fig. [2] we show the variation of the total energy 
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FIG. 2: Variation of total energy (top) and of polarization 
(bottom) as a function of the amplitude of the polar distortion 
between centric (A=0) and polar (A===bl) configurations. 



from the centric to the polar structure as a function 
of the amplitude of the polar mode. For both materi- 
als we find a bi-stable energy profile characteristic of a 
ferroelectric material with an energy barrier on the or- 
der of 0.3eV/unit cell, suggesting that the polarization 
should be switchable upon application of an external elec- 
tric field. For CBDC, the polarization is in the ac plane 
with a mae; nitude of P=14.3/iC/cm^ while for PhMDA 
it is along the c axis and equal to 7.0/iC/cm^. (Recall 
that the polarization of croconic acidis21/iC/cm2). The 
slight discrepancy with experimental values on CBDC 
might be due to insufficient optimization of high-quality- 
crystallization, cite AdvancedMat. We next consider the 
relaxed structures of the high- and low-symmetry phases, 
analyzing the displacive-type transition between the two 
phases in terms of symm etry modes using the Am- 
pHmodes software package |27|. The program deter- 
mines the global structural distortion that relates the 
two phases, enumerates the symmetry modes compatible 
with the symmetry breaking, and decomposes the total 
distortion into amplitudes of these orthonormal symme- 
try modes. 

In Fig. Owe show the centric structures for CBDC (left) 
and PhMDA (right), with the characteristic atomic dis- 
placements of the different polar distortion modes shown 
by colored arrows. As the polar mode acts separately on 
different Wyckoff positions (WPs) of the high-symmetry 
structure, it is meaningful to consider the action of the 
polar distortion on atoms belonging to different WPs sep- 
arately. We denote these as A(WP) and they are shown 
from top to bottom in Fig.O For CBDC, we have A(4e) 
and A{4:b) as intra- and inter-molecular proton transfer 
distortions and A{Sf) as out-of-plane molecular twisting 
(buckling) and pi-bond switching of carboxylic groups. 
For each of them we have calculated the polarization by 
displacing only the atoms belonging to a given WP orbit 
and keeping the rest of them in their centrosymmetric po- 
sitions, obtaining P^]q=(6. 6, 0,-5.5), P4g=(0. 5, 0,-1.6) 
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FIG. 3: (Color online) Displacements patterns (arrows) con- 
necting centric to polar structures for atoms belonging to 
specified Wyckoff positions (top to bottom) for CBDC (left) 
and PhMDA (right). 



and Pg£=(5.5,0,0.4) /iC/cm^. Notably, their sum is 
(12.6,0,-6.7) /iC/cm^ which is almost equal to the total 
polarization P|^q|^ = (12.7,0,— 6.6) /iC/cm^ calculated from 
the total polar distortion. The linearity of the partial 
polarizations is compatible with displacive-type ferro- 
electricity. Surprisingly, the mode decomposition shows 
that the A (8f) mode makes almost as large a contribu- 
tion as the inter-molecular proton transfer in determin- 
ing the total polarization. This effect can be related to 
a double(7r)-bond switching of carboxylic C=0<^=^C— O 
bonds correlated with the "inter-molecular" and "intra- 
molecular" hydrogen distortion. This is shown in Fig. |4l 
where we present an enlarged view of the A(8f) mode. 
In the upper part, we show the centric case, with hydro- 
gens equidistant from nearest carbons or oxygens. In the 
lower part, we show the cooperative hydrogen distortions 
leading to the +P or — P state, which correlates, in turn. 
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do,c, = l-25>^ 

doc=l-25/^ 
dec =1-32 A 



do,c=l-32>^ 

dc,03 = "-25A 
doc = I -32/5^ 
dec =1-25 A 



FIG. 4: (Color online) Switching between double and single 
bonds in the CBDC molecule. Top part: centric case; bottom 
part: switching between +P and —P. The chemist convention 
for the orientation of the dipole moment is used here, i.e. 
arrow starts at 6-\- and ends at 6—. 



with the switching of double and single C-0 bonds and 
with the shortening/elongation of the corresponding C-0 
bonds, in both +P and —P. This picture is further con- 
firmed by the following computational experiment: (i) 
we first consider all the atoms at their centric positions 
(upper part in Fig. S]); (ii) we then move only the intra- 
molecular hydrogen as, for instance, in the +P state, 
keeping all the other atoms fixed. The charge density dif- 
ference between (i) and (ii) (not reported here) shows an 
incipient pile- up of out-of-plane charge between C2=02 
and Ci=03, which corresponds to the incipient forma- 
tion of the TT (double) bonds. One could expect that the 
polar carboxylic groups, rather than the less polar C- 
C bonds, might be responsible for the large polarization 
of the buckling mode. To confirm this, we have further 
extracted the contribution of the C=0/C— O bonds to 
this polarization. The extracted polarization of (4.8, 0, 
-0.1)/iC/cm^ [vs. the remaining contributions of (0.8, 0, 
0.4)/iC/cm^] account for most of the large polarization 
of the A(8f) mode. Again linearity is well fulfilled. This 
clearly explains the origin of the surprisingly large polar- 
ization of the buckling mode. For PhMDA, we find par- 
tial modes 74(4a), A{4:c) and ^4(8^) whose contributions 
to the polarization are 5.8, 1.0 and 0.3 /iC/cm^ along the 
z polar axis, respectively. The linearity holds also in this 
case, but now the inter- molecular proton transfer does 
give the dominant contribution. 

Finally, for the centric structures, we have calculated 
the dimensionless Born (or "dynamical" or "infrared") 
charge tensors Z^^^^ = {Vt/ e)dPa/ duk^^^ where is the 
component of the polarization in direction a and Uk^(3 
is the displacement of atom k in direction /3, 1^ is the 
primitive cell volume, and e is the charge quantum. In 
perovskite ABO3 oxides the ferroelectric tendency is well 



known to be connected with the presence of anomalously 
large Born charges pisi [l9| . It should be noted that in 
low-symmetry cases, as in the present study, the Born 
tensor is not symmetric in its Cartesian indices. There- 
fore, we have split the tensor into symmetric Zg and 
antisymmetric Z\q parts. In the following, we will focus 
on the former, and in particular, on its three eigenvalues 
Ai > A2 > A3 whose physical interpretation is straight- 
forward. Furthermore, only the hydrogen Born tensors 
will be considered. We have also calculated the phonon 
frequencies at the F point; the presence of an imaginary 
frequency usually implies a structural instability, in this 
case of the paraelectric structure. 

Let us first consider CBDC. As expected, the signifi- 
cant deviations of the dynamical tensor with respect to 
the nominal charges involve the "active" H atoms. For 
the "inter-molecular" hydrogen, A^,^* =(3.4,0.4,0.1) and 
for the "intra-molecular" hydrogen, A^^^* =(2.2,0.3,0.3). 
The large values of the Born charges for hydrogens con- 
firm their important contribution to the polarization, as 
found also in other hydrogen-bonded organic molecular 
crystals [28,,^]. The other hydrogens have only negligible 
absolute eigenvalues '^0.1 (some of them becoming neg- 
ative). For the phonons, we found a large non-degenerate 
imaginary F phonon frequency of about 106 cm~^. Ac- 
cording to a symmetry analysis, Q infrared irreducible 
representations exist for all three WP positions with ei- 
ther Au or Bu symmetry. This is not unexpected, as all 
WP positions carry a contribution to the polarization. 
In particular, the eigenvector of the imaginary frequency 
has symmetry B^^ which is polar. After normalization 
to 1 A, we use the Amplimode software for studying the 
corresponding displacement pattern. The largest abso- 
lute where u is the displacement of the atom ac- 
cording to the phonon eigenvector, is 0.29 and 0.35 A 
for "intra-molecular" and " inter- molecular" hydrogens. 
Again, this confirms the dominant role of the two types of 
hydrogen in the ferroelectric properties. For the case of 
PhMDA, we found significant deviations of the dynami- 
cal charge tensor for "inter-molecular" hydrogen, whose 
eigenvalues are Ai^^*=(4.1, 0.4, 0.2). The eigenvalues 
for other hydrogen atoms are smaller than 1. Finally, 
the imaginary phonon frequency is 112 cm~^ with polar 
symmetry B^ Also in this case, the polarization vector 
of the eigenmode has a large displacement of ~ 0.43 A for 
the "intermolecular" hydrogen. 

To summarize, we have studied the origin of ferroelec- 
tricity in CBDC and PhMDA using first-principles calcu- 
lations and symmetry analysis methods. The estimated 
polarizations are as large as ~ 14 /iC/cm^ for CBDC and 
~7/iC/cm^ for PhMDA. We have shown that a partial 
mode analysis is a useful tool for exploring the polariza- 
tion mechanisms. In both compounds, the proton trans- 
fer between (or within) molecular units appears to be the 
main contribution, as confirmed by the large dynamical 
charges and the analysis of the eigenmode displacement 
patterns. However, other contributions, especially the pi- 
bond switching of carboxyhc groups in CBDC may also 
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have a significant weight in the final polarization. Again, 
partial mode analysis has been used to elucidate the ori- 
gin of the unexpectedly large contribution. We hope that 
our study will stimulate further attempts to search for 
new organic ferroelectrics with potentially large electric 
polarizations. 
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